Spin-orbit coupling induced magnetism in the d-density wave phase of La 2 _ a; Ba :r Cu04 
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We study the effects of spin-orbit coupling in the d-density wave (DDW) phase. In the low- 
temperature orthorhombic phase of La2_ a: Ba :E Cu04, we find that spin-orbit coupling induces fer- 
romagnetic moments in the DDW phase, which are polarized along the [110] direction with a con- 
siderable magnitude. This effect does not exist in the superconducting phase. On the other hand, 
if the d-density wave order does not exist at zero field, a magnetic field along the [110] direction 
always induces such a staggered orbital current. We discuss experimental constraints on the DDW 
states in light of our theoretical predictions. 

PACS numbers: 74.20.-z,74.25.Bt,74.20.Mn 



The mechanism of the pseudogap phenomena in high 
T c superconductors remains controversial. Chakravarty 
et al. 1] proposed that it may originate from a hid- 
den long-range (i-density wave (DDW) order 0,13, which 
competes with the c£-wave superconductivity (DSC). This 
scenario has aroused much interest. Extensive analytic 
and numerical investigations have shown its existence un- 
der certain conditions in a variety of one and two dimen- 
sional systems 0, 0, El 13 • However, these states are 
hard to detect experimentally and results are still con- 
troversial 0. Polarized neutron scattering experiments 
0,O] in YBa 2 Cu0 6+x show some supporting evidence. 
On the other hand, Stock et al. ^3 found no indication 
of this phase using non-polarized neutron beams. 

Recently, spin-orbit (SO) coupling has received much 
attention in the emerging science of spintronics. Mu- 
rakami, Nagaosa and Zhang proposed the intrinsic spin 
Hall effect through SO coupling in the p-doped semicon- 
ductors to generate the dissipationless spin current using 
electrical fields ^3 . Similar effects were also predicted in 
the n-doped systems [l4|. The spin-Hall effect in GaAs 
has already been observed experimentally |l5j . 

In Mott-insulators, SO coupling also has important 
effects on the Heisenberg superexchange interactions, 
which is responsible for the anisotropic correction termed 
the Dzyaloshinskii-Moriya (DM) interaction 0]. Given 
the intrinsic spin Hall effect in the semiconductors, it is 
natural to ask what happens in the presence of SO cou- 
pling with states carrying spontaneous electrical currents 
like the DDW state. The answer turns out to depend on 
details, but these conspire in the La 2 Cu0 4 system to give 
rise to an experimentally observable effect: as in the half- 
filled antiferromagnets, SO coupling gives rise to a weak 
planar ferromagnetism which can be used to detect this 
otherwise elusive phenomenon. 



At zero temperature, La2- 2; Ba a: Cu04 undergoes a 
structural phase transition from the low-temperature or- 
thorhombic (LTO) phase to the low temperature tetrag- 
onal (LTT) phase at doping S c s» 0.12. In the LTO phase 
of the undoped La2Cu04, the DM interaction originates 
from the staggered distortion pattern of the oxygen octa- 
hcdrals. It results in the antifcrromagnetic moments ly- 
ing in the afr-plane and the weak ferromagnetic moments 
along the c-axis ^jj- By analogy, the DDW state ex- 
hibits staggered orbital moments. Without SO coupling, 
the DDW state decouples from the spin channel, which 
remains paramagnetic. However, SO coupling couples 
these two channels together, and thus orbital currents 
should affect spin and lead to observable effects. 

In this article, we find that the staggered orbital cur- 
rent induces uniform ferromagnetic moments in both 
LTO and LTT phases. The moments lie in the [110] 
direction in the LTO phase and the [100] direction in the 
LTT phase, respectively. The magnitude per Cu cite is 
at the order of several percent of one Bohr magneton, 
thus is detectable. Conversely, if the DDW order does 
not exist in the ground state, it can be induced by mag- 
netic fields, suggesting that the degree of proximity to 
the instability can be in principle investigated as well. 

We consider the mean field DDW Hamiltonian with 
SO coupling in the LTO phase: 

Hmf = { c L( _i cft + «Ay ■■ a a/ 3)cjp + h. c.j 

(ij)af3 
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where (ij) indicates summation over the nearest neigh- 
bors only. Xij is the decoupling of the Heisenberg ex- 
change term in the particle-hole channel 0] as Xij = 
V{c\ a Cj,j). Its imaginary part Iv&Xij IS the DDW order 
parameter and is treated self-consistently below. On the 
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FIG. 1: A) The lattice symmetry in the LTO phase. 0+ 
(O— ) denotes the oxygen atom moving into (out of) the 
CuO plane. The arrow indicates the [110] direction. B) The 
anisotropic Dirac cone-like dispersion relations around nodes 
(±7r/2, ±7t/2) for spin parallel (|) and anti-parallel (J.) to the 
[110] direction. 

other hand, Rexij changes slowly within the parameter 
regime discussed below, thus is absorbed into the effective 
hopping integral t e f / which gives the band width of holes. 
We choose t e ff sw lOOmeV which is rescaled to 1 below. 
The SO coupling term Ay is determined by the lattice 
symmetries |l8j as shown in Fig. ^ A, such as i) two- 
fold rotations around the c-axes passing the in-plane O 
sites, ii) inversion symmetry with respect to Cu sites, and 
iii) reflection symmetry with respect to the [110] direc- 
tion. Consequently, it shows a staggered pattern Xi^+ X = 
(-)^+MAi,A 2 ,0), X iti+y = (-)^ +l »(-A 2 ,-A 1 ,0) 
Ay is almost perpendicular to the bond direction in the 
LTO phase l2fj|. i. e.. Ai <C A 2 , and A 2 is estimated to be 
around 2meV |l8|. In the LTO phase, the symmetries 
ii) and iii) still ensure the ImXij to exhibit the d x 2_ y 2 
pattern. 

We first present physical arguments for the appear- 
ance of the ferromagnetic moments for a simplified case 
of Ai = A 2 = A/ a/2, and then show that the realistic 
values of Ai j2 essentially give the same result. In the 
simplified case, the SO coupling term plays the role of 
the staggered spin flux with the quantization axis along 
the [110] direction 0]. The effective Hamiltonian reads 

H = J2 C k a C ^k - M) +2i/(fc)4a{ A °X"/3 
k 

+ Imx<W} c fc+Q,/3, (2) 

where f(k) — cosk x — cosk y , o\ = [<j x + a y )/^/2, 
Q = (ir, 7r), and the spin index a =| (1) refers to par- 
allel (anti-parallel) to the [110] direction. For k around 
the nodes (±7t/2, ±7t/2), we define k\\ and k± to be the 
projections of its deviation from the nodes on the di- 
rections parallel and perpendicular to the nested Fermi 
surface, respectively. As shown in Fig. Q]B, we obtain 
the anisotropic Dirac-cone like dispersion relation with 



FIG. 2: Spin polarization in the LTO phase at doping 5 — 0.1 
with Ai = Acos0,A 2 = Asin6> (45° < 6 < 90° ). The results 
within the range 0° < 8 < 45° are symmetric to the case of 
90° - e. 

different slopes for spin f (J.) electrons 

E(k) a , a = ±^v 2 f kl+8(lm X ±\) 2 kl (3) 

where the first ± on the RHS of Eq. @ corresponds 
to the upper (lower) Dirac cone as denoted by the band 
index a, and the second ± corresponds to the spin direc- 
tions | (I), respectively. At half-filling, no ferromagnetic 
moments exist because the lower two bands with opposite 
spin configurations are both fully occupied. However, at 
finite doping S, they are occupied differently, thus a spin 
polarization appears along the [110] direction. At very 
small doping, the ferromagnetic moment per site can be 
estimated from the dispersion relation of Eq. Q as 

M r |A/x| at \X/ X \ « 1 
5 I |x/A] at | X /A| « 1 • {V 

It is instructive to consider the underlying symmetry 
reasons for this effect: the DDW order breaks the time 
reversal (TR) symmetry in the orbital channel while fer- 
romagnetism breaks it in the spin channel. SO coupling 
couples two channels together to linear order. We em- 
phasize that the polarization along the [110] direction is 
valid for general values of Ai j2 . This is protected by the 
following symmetry structures. Although the symmetries 
i) and iii) are broken by the DDW order, their combina- 
tion with the TR operation together still leave the system 
invariant. These symmetries fix the only possible spin po- 
larization along the [110] direction, and further exclude 
the antiferromagnetic order. 

Now we consider the general values of Ai ;2 with the pa- 
rameterization Ai = Acos#, A 2 = Asin^ (9 = 0° ~ 90°). 
The realistic values in La 2 _ 2: Ba 2: Cu04, i.e., Ai <C A 2 , 
corresponds to 6 = 90°. The new effective Hamiltonian 
includes Eq. J5J but with the replacement of A with 
A(cos# + sin 8)/^/2, and also an extra term of 

AH = -2i J2 { AA 9{k)c\ a ^ a0 c u ^ p - h.c.}, (5) 

k,a/3 
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FIG. 3: A) Spin polarizations M v.s. doping levels 8 at 
T=0K. B) The DDW order parameter Imx v.s. doping levels 
5 at T=0K. 
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FIG. 4: Spin polarization v.s. temperatures at S — 0.06,0.12. 
Temperatures are in the unit of t e ff. 
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where g(k) = cosk x + cosk y , er 2 = {cr x — a y )/\/2, and 
AA = A(sin6> — cos#)/\/2. Correspondingly, the spin 
quantization axis of the electron eigenstate depends on 
the momentum fc, which deviates from the [110] direc- 
tion with an angle (j) k satisfying tan0 fc = fr g^fggff ■ 
This helical structure reduces the magnitude of the fer- 
romagnetic moments. However, due to the fact that t e ff 
is much larger than \ an d A, the Dirac-cone is highly 
anisotropic. Then k\\ 3> k±_ holds on most part of the 
Fermi pocket, and the spin deviates from the [110] di- 
rection only at small angles. In other words, the induced 
ferromagnetic moment is insensitive to the ratio of A1/A2. 

In Fig |2 we show the numerical results for the spin 
polarization per site M with the general values of 9 at 
5 = 0.1 and T = OK by using the standard self-consistent 
method [2^]. We choose parameters A = 0.02, V — 0.22 
to agree with the physical value of SO coupling and 
to arrive at a reasonable pseudogap energy scale in 
La2_ a; Ba 2 ;Cu04. In the realistic case of 9 = 90°, the po- 
larization only decreases about 15% compared with its 
maximal value at 9 = 0°. For all numeric results below, 
we keep 9 = 90°, i.e., Ai = 0. We further show M and 
Imx v.s. the doping S at T = OK in Fig. [§Y and|3^. At 
low doping r5, M indeed scales with 6 linearly as indicated 
in Eq. As S increases, Im% drops, and consequently 
M increases faster than linearly. 

The finite temperature behavior of the induced polar- 
ization M is also interesting as shown in Fig. 0] At small 
dopings (e.g. 6 = 0.06), M increases slowly at low T and 



FIG. 5: Staggered orbital moment in the unit of fiB v.s. 
the Zeeman energy E z — g/2B/j,B, where t e // = lOOmeV, 
A/te// = 0.02, AA = 0, and the Lande factor g = 2 . B is 
along the [110] direction. 

decreases after T passes an intermediate value, while M 
decreases monotonically at high dopings (e.g. 8 = 0.12). 
At low dopings, Im% is large at T = OK which decreases 
with increasing T . As a result, the ratio A/Im^ increases, 
and so does M as indicated in Eq. As T goes large, 
then the thermal effect dominates and M decreases. At 
high dopings where Imx is comparable to A, M depends 
less sensitively on A/Imx, and thus thermal effects dom- 
inates in the whole temperature range. 

The magnitude of the induced moments is at the or- 
der of 10 _2 /it_B which should be detectable. Due to their 
ferromagnetic nature, domain wall structures should be 
formed and no macroscopic magnetic field is present. A 
hysteresis behavior appears when an external in-plane 
magnetic field is applied to the sample. In the neutron 
scattering experiment, the elastic Bragg peaks at recip- 
rocal lattice vectors are the evidence for these moments. 
The muon spin relaxation is also sensitive to the internal 
magnetic fields. To our knowledge, no such effects have 
been detected in the LTO phase of I^-ajBazCuCM. 

SO coupling also leads to a "staggered spin galvanic 
effect" as the inverse of the DDW induced ferromag- 
netism. Assuming that the DDW order does not exist 
in the ground state, or equivalently setting V — in 
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Eq. and adding the Zeeman energy term H z = 

— 9^bB - Si, we find that a spin polarization along 
the [110] direction induces a staggered orbital current. 
The staggered current per bond is written as 

ieteff v - ^ + e v— * + 

Istag = -Jfj- 2^ C ka C k+Qa f{k) + — =^ 2^ C k a C W 
ka k,a(3 

X |(Ai COsfca; + A 2 COSk y )a x , a j3 

+ (X 2 cos k x + Xi cos k y )a y ^ a pY (6) 

where the second term originates from SO coupling. Un- 
der the symmetry analysis, I stag can be induced only by 
the B field along the [110] direction. Using the values of 
t e ff and A stated above, the lattice constant a = 3.8 A, 
and = 90°, we show in Fig. El the linear behavior of the 
staggered orbital moment per plaquette v.s. Zeeman en- 
ergy E z — ^gusB. The magnitude reaches the order of 
10~ 3 he at J? 2 ~ lmcV which corresponds to B « 10T. 
The typical value of the DDW orbital moment estimated 
theoretically is at the order of 10 _2 ^b Q. Compared 
with it, our induced orbital moment is about one order 
smaller. 

In the LTT phase, the lattice symmetry results 
in a different staggered SO coupling pattern as 
X^+i = (-)^+ I «(A 1 ,0,0),A M+!} = -H l -+MA 2 ,0,0) 
(19j |. where the spin quantization is fixed along the [100] 
direction and also Ai <C A 2 20]. Similar analysis indi- 
cates that the previous results also apply here with the 
replacement of the [110] with the [100] direction. 

Next we discuss the effect of SO coupling in the su- 
perconducting portion of the phase diagram. SO cou- 
pling does not change the nature of the DSC phase in 
the absence of the DDW order. Due to the conserva- 
tion of the TR symmetry, no ferromagnetic moment can 
appear. Although the spin and the 4-fold rotational sym- 
metries are broken, symmetries ii) and iii) still leave the 
d-wave singlet pairing structure unchanged, and thus the 
node quasiparticles are kept. In the coexisting region 
of DDW and DSC, if the DDW order is large compared 
with the DSC, the induced ferromagnetism suppresses 
the superconductivity, which may lead to Fulde-Ferrell- 
Larkin-Ovchinnikov phases [2lJ . 

At last, we briefly discuss the YBa 2 Cu30 6+( 5 system. 
The inversion symmetry is broken in each CuO plane, 
and the resulting SO coupling is uniform in contrast to 
the staggered pattern in L&2- x Baa;Cu04 systems. Due 
to the existence of CuO chains, the [100] and [010] direc- 
tions are not equivalent any more, and the 4-fold rota- 
tional symmetry is broken, while the reflection symme- 
tries respect to [100] and [010] directions are still kept. A 
straightforward symmetry analysis gives the form of the 
SO coupling H3 as A M+ x = (0, Ai, 0), X lyi+y = (A 2 , 0, 0). 
However, because the staggered orbital current still pre- 
serves the above two reflection symmetries, this SO cou- 
pling can not induce magnetic moments on Cu sites. We 



notice that at least in a single plane the broken parity 
has a significant effect in the superconducting phase: as 
in the 2D Rashba system the singlet and triplet pairing 
channels should be mixed [23j ]. 

In summary, we have investigated the effect of SO cou- 
pling to the DDW state in doped La 2 _ 2; Ca a ;Cu04. If the 
DDW state indeed exists, SO coupling results in a uni- 
form ferromagnetic moment along the [110] direction in 
the LTO phase or the [100] direction in the LTT phase. 
This effect can be used to test the validity of the DDW 
scenario for the pseudogap mechanism. The inverse ef- 
fect is also predicted that a in-plane Zeeman field induces 
a staggered orbit moment. 
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